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' Topological complexity of collision free motion planning 

Qj • algorithms in the presence of multiple moving obstacles 

' Michael Farber, Mark Grant, and Sergey Yuzvinsky 

Abstract. We study motion planning algorithms for collision free control of 
■ multiple objects in the presence of moving obstacles. We compute the topolog- 
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ical complexity of algorithms solving this problem. We apply topological tools 
and use information about cohomology algebras of configuration spaces. The 
results of the paper may potentially be used in systems of automatic traffic 
control. 



1. Introduction 



' The theory of robot motion planning |1U) has developed a broad variety 

I of algorithms designed for different real life situations. Most of these algorithms 

' assume that the robot performs an online search of the scene and eventually finds 

0^ [ its way to the goal. In this approach one understands that the robot has initially 

' a very limited knowledge of the scene and supplements this knowledge by using 

sensors, vision, memory and perhaps some ability to analyze these data. 

A completely different situation arises when one controls simultaneously mul- 
I tiple objects (robots) moving in a coordinated way in a fully known environment. 

■ Assume, for instance, that we have n objects moving in R'^ with no collisions and 
^ I avoiding the obstacles whose geometry is prescribed in advance. In this case the 

■ dimension 3n of the configuration space of the system is large (if the number n of 
. 5^ I controlled objects is large) and therefore the online search algorithms as described 

i above become much less effective. 

}H I The topological approach to the robot motion planning problem initiated in 

[n |5] is applicable in some situations with high dimensional configuration space 
under the assumption that the configuration space of the system is known in ad- 
vance. One divides the whole configuration space into pieces (local domains) and 
prescribes continuous motions (local rules) over each of the local domains. The 
minimal number of such local domains is the measure of topological complexity of 
the problem. 

Formally, the topological complexity is a numerical invariant TC{X) of the 
configuration space of the system [J . Its introduction was inspired by the earlier 
well-known work of S. Smale jl6) and V. Vassiliev 17 on the theory of topological 



1991 Mathematics Subject Classification. Primary 55R80; Secondary 93C83. 
Key words and phrases. Robot motion planning, avoiding obstacles, moving obstacles, 
Schwarz genus, topological complexity, configuration spaces. 

1 



2 



MICHAEL FARBER, MARK GRANT, AND SERGEY YUZVINSKY 



complexity of algorithms of solving polynomial equations. The approach of [Tl l2] 
was also based on the general theory of robot motion planning algorithms described 
in the book of J.-C. Latombe and on the abstract theory of a genus of a fiber 
space developed by A.S. Schwarz |14| . Paper [6] contains a recent survey. 

Paper solves the problem of finding the topological complexity of motion 
planning algorithms for controlling many particles moving in R'^ or on the plane 

with no collisions. Such motion planning algorithms appear in automatic traffic 
control problems. Paper |3] uses the techniques of the theory of subspace arrange- 
ments 12 . It was shown in that the complexity is roughly twice the number 
of particles (it is 2n — 2 for the planar case and 2n — 1 for the spatial case). At 
the moment we do not know specific motion planning algorithms with complexity 
linear in n. However a quadratic in n motion planning algorithm was described in 
§26 of iSj. 

In paper |3] we studied the problem of computing the topological complexity 
of the motion planning problem for n particles moving on a graph with forbidden 
collisions. It was shown that for large numbers of particles the topological complex- 
ity depends only on the geometry of the graph and is independent of the number 
of particles. Some specific motion planning algorithms for this problem were also 
constructed 3.. 

In the present paper we study algorithms solving the following motion plan- 
ning problem. Several objects are to be transported from an initial configuration 
Ai, A2, ■ ■ ■ , An G R"^ to a final configuration Bi, B2, ■ ■ ■ , Bn G R'^ such that in the 
process of motion there occur no collision between the objects and such that the 
objects do not touch the obstacles in the process of motion. 

We make several simplifying assumptions: (a) Each object is represented by a 
single point; (b) each obstacle is represented by a single point; (c) collision between 
two objects occurs if they are situated at the same point in space; (d) an object 
touches an obstacle if the points representing the object and the obstacle coin- 
cide. We also make an important assumption that the behavior of the obstacles is 
known in advance. We show that because of this assumption the problem becomes 
topologically equivalent to the similar problem when the obstacles are stationary. 

One of our results seems surprising: we show that complexity of collision free 
navigation of many objects in the presence of moving obstacles is essentially inde- 
pendent of the number of obstacles and grows linearly with the number of controlled 
objects. 

We believe that the conclusions of this paper will remain valid in a more general 
and realistic situation when the objects and the obstacles are represented by small 
balls, possibly of different radii, and the control requirements are to avoid tangencies 
between objects and obstacles. 

2. Statement of the result 

Let us explain our notation, n denotes the number of moving objects and m 
the number of obstacles. The symbol F{X^ n) denotes the configuration space of n 
distinct points in a topological space X; in other words, F{X, n) is the subset of the 
Cartesian power X x X x ■ ■ ■ x X (n times) consisting of configurations {xi, . . . , a;„) 
(where Xi G X) with Xi 7^ Xj for i ^ j- 

We assume that the controllable objects lie in R*""*"^ where r > 1 is an integer. 
Clearly, the practical applications are covered by the cases r = 1 or ?■ = 2 which 
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will be considered in more detail later. In this section we will assume that r is 
arbitrary. 

The algorithms we study take as input the following data: 

(a) the initial configuration A — {Ai, A2, . . . , A^) £ F(W^^,n) of the ob- 
jects, Ai e W+^, with Ai ^ Aj] 

(b) the desired final configuration of the objects B = {Bi, B2, ■ ■ ■ , Bn) G 
F{W+\n), where Bi G B, ^ By, 

(c) the trajectory of moving obstacles C{t) — {Ci{t),C2{t), . . . ,C„i{t)) G 
F(R''+i,m), Q{t) ^ Cj{t) for any t G [0, 1]. 

Here t G [0,1] is the time. We assume that Ai / C'j(O) and Bi ^ Gj{l) for all 
i = 1, . . . , n and j = 1, . . . , to. This means that initially the objects are all distinct 
and disjoint from the obstacles at time t = and the desired positions of the objects 
are pairwisc distinct and distinct from positions of the obstacles at time t = 1. 
The output of the algorithm is a continuous motion of the objects 

-f,{t)eIC+\ iG[0,l], 1=1,..., n 

satisfying the following conditions: 

(a) 7i(i) are continuous functions of t; 

(/?) 7^(0) = A, and 7,(1) = B,; 

(7) for any t G [0, 1] one has 7i(t) ^ Jkit) for i ^ A: where z, fc G {1, . . . , n}; 

(6) for any t G [0, 1] one has 7i(i) 7^ Cj{t) where i = 1, . . . ,n and j = 1, . . . ,m. 

Let the function 

: [0,1] ^F(R'^+\m) 

represent the motion of the obstacles. For any time t G [0, 1] the obstacles are 
positioned at the points 

Y^{t)^{Ci{t),C2{t),...,an{t)). 

Next we introduce a few more notations: 

E{Ym,) = {7 : [0, 1] ^ F{W+\ny, j{t) n Y„,{t) =0, Vt G [0, 1]}, 

B{Y,n) - {{A,B) G F{W+\n) x F(R'-+\ n); ^ n r,„(0) - 0, BnY„,{l) = 0}. 

Elements of E{Ym) are motions of the objects avoiding the obstacles, elements 
of B{Ym) represent the input data: the initial configuration A and the desired 
configuration B. 

There is a canonical map 

(2.1) ^(y„) : E{Y^) -> B(y„) 

given by 7r(y„i)(7) — (7(0), 7(1)). A section of tt is a map (possibly discontinuous) 

(2.2) s : B{Y^) ^ E{Yrn), 

satisfying Tr(Ym) o s = lg^Y„^)- Clearly, s assigns to any input data (A, B) G B{Ym) 
a motion of the controllable objects starting at configuration A, ending at B and 
avoiding mutual collisions between objects and between objects and the moving 
obstacles. We conclude that algorithms solving our motion planning problem are 
in one-to-one correspondence with sections of map (j2.1|l (we will see below that it 
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is a fibration). Given a section s as above, then for any pair {A,B) E B{Ym), the 
output of the algorithm is s{A,B) E E{Ym)- 

Hence, our task is to estimate the complexity of finding a section of iT{Ym)- 
We recall that the Schwarz genus 14 of a continuous map f : X ^ Y map is 
the minimal integer k such that Y can be covered by k open sets UiL)U2 ■ ■ -UUk ^ Y 
with the property that over each Ui there exists a continuous section Si : Ui —> X, 
f o Si = Ijji- We apply this notion to estimate the complexity of finding a section 
to l(TT|l . 

Definition 2.1. The complexity of the motion planning problem for moving 
n objects in the presence of m moving obstacles is defined as the Schwarz genus of 
map (|2.1|) . 

In practical terms this means the following. One divides the space of all pos- 
sible inputs B{Y„i) into several domains (called the local domains) and specifies 
a continuous motion planning algorithm (called a local rule) over each of the do- 
mains. The complexity of the problem is defined (according to Definition 12.1(1 as 
the minimal number of local domains in all possible algorithms of this type. 

Now we state our main result: 

Theorem 2.2. The complexity of motion planning problem of moving n > 2 
objects in avoiding collisions in the presence of m > 1 moving obstacles equals 
2n + 1. The complexity of the similar planar problem (i.e. when the objects are 
restricted to lie in the plane R^) is 2n if m = \ and it is 2n+ 1 for m > 2. 

This theorem follows by combining Theorem 13.11 and Theorems 15.11 and 16.11 
which are stated and proven below. 

3. Reduction to the case of stationary obstacles 

The following theorem shows that Definition 12.11 gives the notion of complexity 
which coincides with a special case of the notion of navigational complexity of 
topological spaces TC{X) which was studied previously [Jj. Recall that TC{X) 
denotes the Schwarz genus of the path space fibration 

(3.1) tt:X^->XxX 

where X^ denotes the space of all continuous paths j : I — [0,1] ^ X with the 
compact-open topology and n{-f) — (7(0), 7(1)). 

Theorem 3.1. (a) The map T:{Ym) (given by 1^2. ) is a locally trivial fibration; 
(b) The fiberwise homeomorphism type of fibration i^iYm) is independent of the 
trajectory of moving obstacles Ym', in particular, it equals the Schwarz genus of the 
special case of 1^2.1]) when the obstacles are stationary; (c) The Schwarz genus of 
T^{Ym) equals 

TC(i^(R''+i-5™,n)) 
where Sm — {si, • • ■ , Sm} C R*""*"^ is an arbitrary fixed m-element subset. 

Proof. First we find a continuous family of homeomorphisms 

V-t : R"+' ^ R''+\ te[0,i] 

having the properties 

(3.2) ^0 = id, and MYmit)) = YniO), t E [0, 1]. 
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Equation (|3.2|l is understood as an equality between ordered sets; in other words, 
(|3.2(l states that 

(3.3) MC^{t)) = Cm. ie[0,l]. 

The existence of such isotopy "0* is a standard fact of the manifold topology; formally 
it follows from the well known isotopy extension theorem. 

Next we consider the commutative diagram of continuous maps 

E{Y,n) ^ X' 



(3.4) i 7r(y„) i TT 

B{Y^) ^ XxX 

where 

X = F(R'-+i-y,„(0),n) 

and TT is the canonical path space fibration H3.1|l for the special case X — F(R''"+^ — 
Y„i{Q),n). Maps F and G are defined as follows: 

(3.5) F{j){t) ^ Mlit)), l^E{Ym), te[0,l], 



(3.6) G{A,B) = {MA),MB)), {A, B) e B{Y„,). 

Clearly, F and G are homeomorphisms and diagram (|3.4(l is commutative. This 
implies statements (a) - (c) of Theorem 13. II □ 



4. Cohomology of configuration space i^(R''+^ — Sm,n) 

In this section we collect some known topological results which will be used in 
this paper. All cohomology groups have Z as coefficients. 

Let Sm = {si, . . . ,Sm} C R*""*"^ be a fixed sequence of m distinct points in 
Euclidean space R''+^. Here r > 1 is an integer. The space F(R''+^ — Sm, n) repre- 
sents configurations of n particles in R''+^ avoiding mutual collisions and collisions 
with m stationary obstacles si, . . . , Sm- 

The theorems we state in this section are known, they can be found in Chapter 
V of the book by Fadell and Husseini [3; we have adjusted the notations as required 
for our needs. The space i^(R^ — Sm,n) (here r = 1) can be identified with the 
complement of an arrangement of affine hyperplanes in C", so one can also use the 
book | 12| as the reference. 

We begin by noting that F(R''+^ — Sm,n) embeds in the configuration space 
F{W+'^,n + m) via the map (?/i, . . . ,y„) i-^ (yi, . . . , y„, si, . . . ,s„i) (in fact this em- 
beds F(R''+^ — Sm, n) as the fibre over (si, . . . , Sm) of the locally trivial fibration 
FiTV^^ -I- m) — > F(R''+^, m) which projects to the last m points of the config- 
uration). The cohomology algebra 77*(F(R''+^, n -f m); Z) is well known, and can 
be described as follows. 

For each pair of integers i ^ j such that 1 < i, j < n + m, consider the map 

4j :F(R'■+^^ + ^)^^^ (yi,y2,...,y„+™)^r^^ ^ s\ 

\y^-y3\ 
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Letting [5""] G H'''{S^) denote the fundamental cohomology class of the sphere 5"", 
we obtain cohomology classes 

e,j- = (/)*j-[S'n e H''{F{W+\n + m); Z). 

Theorem 4.1. (7 , Theorem V.4-2) H* {F{IV'^'^ ,n+m)) is the free associative 
graded- commutative algebra generated by the classes Cij for 1 < i < j < n + m, 
subject to the relations 
(i) ef^ = 0, and 

(a) e^jCik = (e.y - eik)ejk for any triple i < j < k. 

Theorem 4.2. (7^, V.4-2) The homomorphism 

H*{F{W+^,n + m)-Z) ff*(F(R'^+i - S,„,n);Z) 

induced by inclusion is an epimorphism, with kernel equal to the ideal generated by 
those Cij having i > n and j > n. 

Corollary 4.3. In H*{F{IU'^'^ - Sm,n); Z) there are relations 

for any triples i,j, k such that i < n and j, k > n. 

Theorem 4.4. (7 , Theorem V.4.3) A basts for H*{F{W+^ - S„„n);Z) is 
given by the set of monomials 

Biiji S-i2j2 ■ ■ ■ ^isjs 

where 1 < ii < Z2 < • ■ • < is < ?i and for each q with 1 < q < s we have iq < jq 
and 1 < jq < n + m. 

We remark here that the maximum length of such a monomial is (n — 1) when 
m = 0, and n when m > \. Hence the highest dimension in which the integral 
cohomology of FiTV'^^ — Sm, n) is non-trivial is r{n — 1) when to = and rn when 
TO > 1. Note also that when m > 1 the given basis includes all monomials of the 
form 

(4.1) /^-f = n "^^^'i+i n '^«.«+2 

for every J C {1, . . . , n}. For r — 1 this fact also follows from the notion of nbc- 
basis for an appropriate ordering of the generators (see \V2\ ). 

In what follows we will specialise to the cases r = 2 and r = 1. 

5. Computing TC(F(R3 - Sm,n)) 
Theorem 5.1. One has: 

2n — 1 if m = 0, 



TC(F(R^-5,„,n)) 



2ri + 1 if m>l. 



Proof. We first establish a lower bound. We shall use the cohomological lower 
bound given by Theorem 7 of jlj . Set = 1 ® Cij — Cij (X) 1 . It is a zero-divisor of the 
cohomology algebra. Note that (cij)^ — —2{eij eij) ^ 0. Consider the following 
product TT = \^Z\ ^inf ■ We find tt = (— 2)"~^jU ® ji, where fi = n"=/ ^in- The 
monomial /z is nonzero by Theorem l4.4l and hence the product tt of length 2(n — 1) 
is nonzero. This gives the lower bound TC(i^(R'^ — S',„, n)) > 2n — 1. Now assuming 
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that TO > 1 we have a nontrivial product Y[7=i of length 2n, which gives 

the lower bound TC(F(R3 - Sm,n)) > 2n + I when to > 1. 

To obtain the upper bound, note that i^(R'^ — Sm,n) can be viewed as the 
complement of a finite collection of codimension 3 affine subspaces in R'^", so it is 
simply-connected by an easy transversality argument. Since it has finitely generated 
torsion-free homology and cohomology, it has the homotopy type of a CW-complex 
consisting of one fc-cell for each fc-dimensional element in the basis for cohomology 
given by Theorem l4.4l fsee Chapter 4.C of [S]). This minimal cell structure is made 
explicit in Theorem VI. 8. 2 of In particular, it has the homotopy type of a 
polyhedron of dimension 2(n — 1) when to = and 2n when m > 1. We now apply 
Corollary 5.3 of 2^ stating that for a 1-connected polyhedron Y, 

TC{Y) < dim(r) + 1, 

which together with homotopy invariancc of TC completes the proof. 

□ 

6. Computing TC(i^(R2 - S,n,n)) 
Theorem 6.1. One has: 

{2n-2 i/ m = 0, 
2n if 1, 

2n + 1 if m>2. 

Proof. The statement for the first two cases follows immediately from |5] and 
the fact that F(R^ — Xi, n) is homotopy equivalent to F(R^, n+1) (see [Z1,P-15). 

In what follows we assume that to > 1. We first establish a lower bound using 
again the cohomological lower bound given by Theorem 7 from [l]. Set as above 
Cij = 1 ® Cij — eij ® 1 and consider the product 

i— l,...,n, J— n+l,n+2 

It is clear that tt can be expressed as a linear combination of pure tensors /ii ® /i2 
where /i^ are monomials in the Cij complementary to each other in tt. Since the 
highest non-zero dimension in H*{F(R?~Sm, n); Z) is n both /ii should have degree 
n in order for /ii ® /i2 not to vanish. Also it follows from the relation of Corollary 
14.31 that the nonvanishing summands of tt are of the form fij (E> fif, where fii is the 
monomial defined by (|4.1|l . / runs over all subsets of {1, 2, . . . , n} and / denotes 
the complement of /. Since the set {fJ,i\I C {!,..• ,n}} is a subset of a basis of 
i7"(i^(R^ — Sm, n); Z) given by Theorcm l4.4l no cancellations are possible, whence 
TT^O. This gives the inequality TC{F{Il'^ - 5„, n)) >2n + l. 

The opposite inequality follows immediately from Theorem 5.2 of 2 since 
F(R^ — Sm^n) has the homotopy type of a connected polyhedron of dimension n 
(see |12p. This completes the proof of the Theorem. □ 

7. Concluding remarks 

First we note that our main result Theorem 12.21 follows by simply combining 
Theorems O O and lO 

Let us compare the following two control problems: (1) motion planning for 
moving n objects in R'^ with no collisions and avoiding collisions with m > 1 
moving obstacles; and (2) motion planning for moving n objects in R'^ with the 
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only requirement that they avoid a single point obstacle (in particular, the objects 
are allowed to collide, i.e. to occupy the same position in space). Note that the 
integer m in problem (1) can be arbitrarily large. According to Theorem 12.21 the 
topological complexity of the problem (1) is 2n+ 1. Surprisingly, problem (2) also 
has complexity 2n+l. Indeed, the configuration space of problem (2) is Cartesian 
power of n copies of — {0} which is homotopy equivalent to 

5*2 X X •■• X (n times) 

which has topological complexity 2n + 1 as it is easy to see. Hence, surprisingly 
problem (1) which is intuitively more "complicated" has the same topological com- 
plexity as problem (2). 

This comparison shows that in general the notion of topological complexity is 
only a partial reflection of real difficulty of a motion planning problem. 
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